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Beschreibung 

Verfahren zur Bilderstellung fur einen mehr zeiligen Spiral- 
Computertomographen und Computer-Tomographie-Gerat zur Durch- 
ftihrung dieses Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bilderstellung fur 
einen mehrzeiligen Spiral - Comput ertomographen , wobei ein Un- 
tersuchungsobjekt durch ein rotierendes und sich in Richtung 
der Rotationsachse bewegendes Strahlenbundel beziiglich seines 
Absorptionsverhaltens abgetastet wird und die gemessenen Ab- 
sorptionsdaten gesammelt werden. Zur Rekonstrukt ion eines Vo- 
lumenbildes aus den gemessenen Daten werden diese gefiltert 
und die gefilterten Daten anschlieftend zur Erzeugung eines 
Volumenbildes des Untersuchungsob j ekts dreidimensional ruck- 
projiziert, wobei das Volumenbild aus den Daten gewonnene Ab- 
sorptipnswerte der zu dem Volumen des Untersuchungsob j ekts 
gehorigen Voxel fur die Strahlung des Strahlenbundels repra- 
sentiert. Aufterdem betrifft die Erfindung ein CT-Gerat zur 
Durchfiihrung dieses Verfahrens mit ent sprechenden Mitteln zur 
Bilderstellung . 

Es ist bekannt, in der Computertomographie Mehrzeilendetekto- 
ren und Strahlenbundel mit Kegelstrahl-Geometrie zur Abtas- 
tung von Untersuchungsobj ekten, insbesondere von Patienten, 
zu verwenden. Zur Rekonstruktion eines Volumenbildes, beste- 
hend aus einer Vielzahl von kleinen Volumenelementen 
(=Voxel) , muss die Kegelstrahl -Geometrie in der nun dreidi- 
mensional en Bildrekonstruktion berticksichtigt werden . 

Solche Rekonstruktionsalgorithmen zur Mehrzeilen- Spiral -CT 
konnen in die zwei Klassen der approximativen Algorithmen und 
der exakten Verfahren gegliedert werden. Beziiglich der appro- 
ximativen Algorithmen wird auf die Schriften : 

- M . Kachelrieii, S. Schaller, and W. A. Kalender, "Ad- 
vanced single-slice rebinning in cone-beam spiral CT" , 
Med. Phys. 27 (2000) 754-772 
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- S. Schaller, K. Stierstorfer, H. Bruder, M. Kachelriefi, 
and T. Flohr, "Novel approximate approach for high- 
quality image reconstruction in helical cone beam CT at 
arbitrary pitch", Proceedings SPIE 4322 (2001) 113-127 
5 - K. Stierstorf er , T. Flohr, H. Bruder, "Segmented Multi- 

ple Plane Reconstruction - A Novel Approximate Recon- 
struction Scheme for Multislice Spiral CT" , Proceedings 
of Intern. Meeting on Fully 3-D Image Reconstruction in 
Radiology and Nuclear Medicine, Pacific Grove, CA, USA, 
10 10/30-11/2/2001, pp. 95-97 

verwiesen. Eine Ubersicht der exakten Verfahren ist in der 
Schrift K. Sourbelle, H. Kudo, G. Lauritsch, K. C. Tarn, M. 
Defrise, and F. Noo, "Performance Evaluation of Exact Cone- 
Beam Algorithms for the Long-Object Problem in Spiral Com- 
- 15 puteci Tomography" , Proceedings of Intern. Meeting on Fully 3- 
D Image Reconstruction in Radiology and Nuclear Medicinev Pa- 
cific Grove, CA, USA, 10/30-11/2/2001, pp. 153-156, aufge- 
fuhrt . 

20 Approximative Verfahren zeichnen sich durch ein hohes Mafl an 
Praktikabilitat und Flexibilitat aus . Allerdings wird der 
Neigungswinkel der Messstrahlen zur Rotationsachse (Kegelwin- 
kel) nur naherungsweise ■ berucksichtigt , weshalb der Approxi- 
^ mationsfehler mit dem Kegelwinkel wachst. Insgesamt lasst 

is sich sagen, dass j edes approximative Verfahren ab einer ge- 



wissen Anzahl von Detektorzeilen Bildartef akte verursachen 
wird und damit zu unbef riedigenden Bildergebnissen flihrt. 
Die exakten Verfahren berucksicht igen die kegelstrahlartige 
Auf nahmegeometrie fehlerfrei sowohl im Filterschritt als auch 
30 in der 3D Riickpro j ektion . Sie erzielen gute Bildergebnisse, 
die unabhangig von dem auftretenden Kegelwinkel sind. 
Nachteilig ist jedoch, dass sie extrem aufwendig und sehr un- 
flexibel in der Anwendung sind. 



35 



Die Aufgabe der Erfindung liegt darin, ein Verfahren zur Bil- 
derstellung fur einen mehrzeiligen Spiral-Computertomographen 
zu finden, welches einerseits die guten Eigenschaf ten der ex- 
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akten Verfahren bezuglich der Bildqualitat beibehalt, ande- 
rerseits jedoch durch Naherungen den Berechnungsauf wand zu 
verringern und moglichst auch die Flexibilitat der Berechnung 
zu erhohen. 

Diese Aufgabe wird durch die Merkmale der unabhangigen Pa- 
tentanspriiche gelost. Vorteilhafte Weiterbildungen der Erf in- 
dung sind Gegenstand untergeordneter Anspruche . 

Exakte Verfahren konnen als Algorithmus vom Typ gefilterte 
Ruckpro jektion (FBP) formuliert werden. Beispiele hierfur 
sind in den Schriften: 

- H. Kudo, F. Noo, and M. Defrise, "Quasi -Exact Filtered 
Backproj ection Al gori thm for Long -Ob j ec t Probl em in Hel i - 
cal Cone-Beam Tomography" , IEEE Transactions on Medical 
Imaging 19 (2000) 902-921; 

- K.C. Tarn, "Exact local regions -of -interest reconstruction 
in spiral cone-beam filtered-backproj ection CT: theory", 
Proceedings of the SPIE 3979 (2000) , 506-519; 

- G. Lauritsch, K.C. Tarn, K. Sourbelle, and S. Schaller , 
"Exact local regions -of -interest reconstruction in spiral 
cone-beam filtered-backproj ection CT: numerical implemen- 
tation and first image results" , Proceedings of the SPIE 
3979 (2000) , 520-532; 

dargestellt . 

Bei diesen exakten Verfahren wird ein Detektorbild X durch 
eine Filteroperation in ein gefiltertes Detektorbild Y iiber- 
fiihrt und anschlieftend mittels 3D Ruckpro j ekt ion das gefil- 
terte Detektorbild Y in ein Volumenbild ubertragen . 

Im folgenden wird ohne Beschrankung der Allgemeinheit der 
Fall eines planaren Detektors beschrieben, der von einem Ko- 
ordinatensystem (u,v) aufgespannt wird. Die v-Achse sei pa- 
rallel zur Drehachse einer Spiralbahn, auf der sich die 
Strahlenquelle und der Detektor bewegen, und die u-Achse 
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■0. 



15 



20 



senkrecht zur dieser Drehachse . Der Filterschritt besteht 
dann aus 

rj+ff/2 

Y=-D ? $d0B(0)P(0)MD HO) X + Bt (1) 
mit 

D- als partielle Ableitung in Richtung des Vektors 

e , 

t als Projektion der Spiraltangente auf den Detek- 

tor, 

als Winkel des Vektors t mit der u-Achse des De- 
tektorkoordinatensystems , 

als Projektion und Ruckpro j ektion auf dem Detek- 
tor entlang Geraden mit 
als Normalenvektor , 

als der Winkel des Normalenvektors r(0) mit der 
u-Achse des Detektorkoordinatensystems , 
als Maskierungsoperation und 

als ein Paar von Maskierungslinien, so dass. in 
der Maskierungsoperation Detektorpixel , die nicht 
zwischen den Maskierungslinien V liegen, zu Null 
gesetzt werden, und 
Bt als Grenzterm (boundary term) . 



P(9) , B(9) 
r{0) 

e 

M 
V 



' '25 
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Die Maskierungslinien V ergeben sich aus der Kegelstrahlpro- 
j ektion der Spiralbahnsegmente , die einen Vollumlauf vor und 
nach der aktuellen Fokusposition liegen. Die Maskierungslinie 
Voben aus der Projektion des Spiralsegmentes oberhalb der ak- 
tuellen Fokusposition schreibt sich in der Form der v- 
Detektorkoordinate als Funktion der u-Detektorkoordinate als 



Y oben ( M ) = "I k 

2p 



V 5 J 



f — - — arctan — | 
U n 5) 



(2) 



mit 
5 

P 

h 



als Abstand Fokus zu Detektor, 
dem Fokusbahnradius , 

dem Liegenvorschub pro Spiralumdrehung . 



\ 
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Die Maskierungslinie Vunten aus der Projektion des Spiralseg- 
mentes unterhalb der aktuellen Fokusposition ergibt sich aus 
der Punktspiegelung von V 0 ben am Detektorkoordinatenursprung 
V unt M = ~V oben {-u) . (3) 

5 Der Grenzterm Bt ist gegeben durch 

rj+x/2 -ko 

Bt = D- \ d6 B(0) \ds(0) (D IW M)s{9) X ( 4 ) 

rj-n/2 -co 

mit dem Vektor s(0) senkrecht zu dem Normal envekt or r(0) der 
Proj ektionsgeraden, s (0) der Koordinatenwert des gedrehten 
Detektorkoordinatensystems (r,s) (0) aufgespannt von den Vek- 
10 toren r(0) und s(0) . 



15 



20 



Unter Ausnutzung der Produktregel flir Ableitungen 
D H0) (MX) = MD fm (x)+D H0) (M)X (5) 

urid der Identitat 

R 7 =-D 7 \dOB{0)P(O)D He) (6) 

TJ-ft/2 

mit R 7 als Rampenfilter in Richtung t der Projektion der 

Spiraltangente wird die Gleichung (1) zu 

Y = R 7 MX + Bt + C (7) 

mit Korrekturterm C als 

rj+x/2 

C = D 7 \d0B(0)P(0)(D He) M)x . (8) 

rj-x/2 



w?* Betrachtet man die Terme der Gleichung (7) , so ergibt sich 
folgendes Bild: 

• Der Term R 7 M X ist schnell und einfach mit Hilfe bekann- 

25 ter Faltungsverf ahren zu berechnen. Er ist in axialer 

Richtung lokal begrenzt . 

• Der Grenzterm Bt ist einfach zu berechnen. Er hat in axia- 
ler Richtung unendliche Ausdehnung, so dass das gefilterte 
Detektorbild Y in solch einer axialen Ausdehnung berechnet 

3 0 werden muss, dass die Kegelstrahlproj ektion des Volumen- 

bildes vollstandig darin enthalten ist. Dies fuhrt zu ei- 
nem erhohten Aufwand in der 3D-Ruckproj ektion . Zudem wer- 
den dadurch Objektpunkte des Volumenbildes kunstlich kor- 
reliert, was zu dem bekannten "Long Object Problem" fiihrt. 



\ 
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• Der Korrekturterm C ist numerisch sensitiv und somit auf- 
wendig zu berechnen. Hier wird ein Oversampling benotigt, 
was die Rechenzeit in die Hohe treibt. Der Korrekturterm 
hat in axialer Richtung unendliche Ausdehnung. Die Konse- 
5 quenzen wurden bereits oben beim Grenzterm Bt diskutiert. 

Die Erfinder haben nun erkannt , dass es zu einer wesentlichen 
Verbesserung der Praktikabilitat des oben beschriebenen exak- 
ten Verfahrens fiihrt, wenn der Grenzterm Bt und der Korrek- 
10 turterm C vereinfacht werden kann. Insbesondere kann in axia- 
ler Richtung die unendliche Ausdehnung vermieden werden. 

Durch einen glucklichen Griff der Erfinder hat sich gezeigt, 
dass der Term 2* (Bt+C) naherungsweise einem Term entspricht, 
15 der wechselweise die Maskierung M und Filterung R 7 auf das 

Detektorbild X anwendet und anschlieiiend die Differenz bil- 
det . Es gilt also : 

Bt + C* -(MR r - R 7 M)X . (9) 
2 

Setzt man Gleichung (9) in Gleichung (7) ein, so erhalt man 
2 0 die einf ache Filtervorschrif t 

y = ^{r 7 m+mr 7 )x . (10) 



Durch die Anwendung dieser Filtervorschrif t lassen sich die 



M folgenden Vorteile erreichen: 



TP 

2 5 • Die beiden Terme von Gleichung (10) sind einfach und 

schnell mittels Faltungsverf ahren zu berechnen. 

• Das gefilterte Detektorbild Y ist in axialer Richtung lo- 
kal begrenzt . Neben den Rechenzeiteinsparungen in der 3D 
Ruckproj ektion umgeht man damit auch das "Long Object 

3 0 Problem" . 

• Die Approximation liefert bis zu relativ hohen Kegelwin- 
keln, entsprechend einer grofien Anzahl an Detektor zeilen f 
gute Bildqualitat . 

• Die Ausnutzung von Datenredundanzen bei reduzierten, dis- 
35 kreten Werten des Liegenvorschubs ist wesentlich einfacher 

als bei exakten Algorithmen. Die Approximation macht , wie 



\ 
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spater erlautert wird, eine unterschiedliche Gewichtung 
mancher Beitrage unnotig. 

Zur Durchfuhrung des Verfahrens sollte beachtet werden, dass 
5 die Anwendung des Rampenf ilters R r auf ein maskiertes Detek- 

torbild MX eine Glattung des Ubergangs vom Innenbereich der 
Maske mit Wert 1 zum AufJenbereich mit Wert 0 verlangt, da der 
Rampenf ilter an unstetigen Stellen Bildartef akte generiert . 
Zur Vermeidung dieser Bildartef akte sollte die geglattete 
10 Maske in Richtung t des Rampenf ilters zumindest weder Unste- 
tigkeiten in der Maske selbst, noch in der ersten, partiellen 
Ableitung aufweisen. Vorteilhaft kann fiir einen glatten Uber- 
gang eine sin 2 -Funktion verwendet werden, wobei darauf hinzu- 
weisen ist, dass auch andere ahnliche Glattungsf unktionen 
15 verwendet werden konnen. 



In medizinischen Anwendungen ist es manchmal vorteilhaft, den 
Liegenvorschub zu reduzieren. Aus der Messung und Verrechnung 
redundanter Daten soli bei reduziertem Liegenvorschub das 
20 Signal-zu-Rauschen Verhaltnis vergrofiert werden. Weder der 

exakte Algorithmus nach Gleichung (7) noch die Naherung nach 
Gleichung (10) bieten unmittelbar die Moglichkeit der Ausnut- 
zung von Datenredundanzen, da die genutzte Detektorf lache bei 
reduziertem Liegenvorschub gemaii der Maskierung nach Glei- 
chung (2) und (3) zusammenschrumpf en wurde . 

Fur bestimmte diskrete Werte des Liegenvorschubs konnen bei 
exakten Berechnungsverf ahren Datenredundanzen genutzt werden. 
Wird der maximale Liegenvorschub h max als der Liegenvorschub 
3 0 definiert, bei dem die durch die Maskierung M ausgezeichnete 
Flache noch von dem Detektor abgedeckt wird, so konnen bei 
diskreten Liegenvorschubwerten hk, 

h. = hmax , mit natiirlicher Zahl k > 0 (11) 
* 2k + 1 

die vorhandenen Datenredundanzen genutzt werden. Dazu wird 
35 der maskierte Bereich erweitert . Die Maske M k wird definiert 
durch die Kegelstrahlproj ektion der Spiralbahnsegmente der 
(k+l)-ten Umlaufe vor und nach der aktuellen Fokusposition . 



\ 
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1 1 ♦ «0 
— - — arctan — . 

In 8) 



V kjob M = — -K • I+-5- 1-| * + --- arctan T | . (12) 



Die Maskierungslinie V k(0ben des Spiralbahnsegmentes oberhalb 
der aktuellen Fokusposition ist gegeben durch 

S u fi " 2 
h. • 1 + —, , , 

Fur die Maskierungslinie V k/Unten gilt die Symmetrierelat ion 
5 von Gleichung (3) . Gleichung (11) in Gleichung (12) ergibt 

Knt>, n (u) = * +-t N arctan— . 13) 

Die Maske M k ist immer vollstandig von jener rechteckigen De- 
tektorf lache abgedeckt, die auch die Maske M=M 0 des maximalen 
Liegenvorschubes vollstandig • beinhaltet . Betrachtet man das 
jj^p Volumen, das von der 3D-Ruckprojektion solcherart maskierter 
^ Detektorbilder uberdeckt wird, so sieht man, dass das Volu- 

menbild naherungsweise (2k+l)-mal uberdeckt wird. Die Filter- 
vorschrift von Gleichung (10) wird auf einfache Weise gean- 
dert zu 

15 Y =- — (R 7 M k +M k R 7 )x . (14) 

Entsprechend dem oben ausfiihrlich dargestellten Grundgedanken 
der Erf indung schlagen die Erf inder vor, das an sich bekannte 
Verfahren zur Bilderstellung fur einen mehrzeiligen Spiral- 

2 0 Computertomographen zu verbessern. Bekannt ist das Verfahren, 
soweit ein Untersuchungsob j ekt durch ein rotierendes und sich 
in Richtung der Rotationsachse bewegendes Strahlenbiindel be- 
^ ziiglich seines Absorptionsverhaltens abgetastet wird und die 
gemessenen Absorptionsdaten gesammelt werden, wobei zur Re- 

2 5 konstruktion eines Volumenbildes aus den gemessenen Daten 
diese gefiltert und die gefilterten Daten anschliefiend zur 
Erzeugung eines Volumenbildes des Untersuchungsob j ekts drei- 
dimensional ruckpro j iziert werden, wobei das Volumenbild aus 
den Daten gewonnene Absorptionswerte der zu dem Volumen des 

30 Untersuchungsobj ekts gehbrigen Voxel fur die Strahlung des 

Strahlenbundels reprasentiert . Die erf indungsgemalie Verbesse- 

rung liegt nun darin, dass die Filterung durch mehrfache An- 
wendung eines Rampenf ilters R 7 und einer Maskierungsoperati- 
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on M auf ein Proj ektionsbild in unterschiedlicher Reihenfol- 
ge erf olgt . 

Gemali einer vorteilhaf ten Ausbildung des Verfahrens wird das 
Rampenfilter R 7 und die Maskierungsoperation M abwechselnd 

angewendet, wobei bevorzugt das gefilterte Detektorbild Y mit 
der Gleichung 

Y =~(R 7 M + MR 7 )X 

berechnet wird, worin 

X das ungefilterte Detektorbild, 

R 7 die Rampenf ilterung in Richtung der Proj ekt ion der Spi- 

ral tangent e und 
M die Maskierung darstellen. 

Gemaft einer weiteren besonderen Ausflihrung des Verfahrens 
schlagen die Erfinder vor, dass ein Liegenvorschub h k fur die 
Relativbewegung von Fokus und Untersuchungsobj ekt in Richtung 
der Rotationsachse des Fokus verwendet wird, dessen Wert der 
Bedingung 



geniigt, wobei k eine natiirliche Zahl mit k > 0 ist und h max 
dem Liegenvorschub entspricht, bei dem eine durch die Maskie- 
rung ausgezeichnete Flache gerade noch den Detektor abdeckt. 
Durch diese Bedingung kann bei einem linearen Liegenvorschub 
mit den diskreten Werten h k die Filterung mit der Vorschrift: 

Y = — (R 7 M. + M.R 7 )X 

durchgefuhrt werden, wobei k die Indexzahl aus dem Liegenvor- 
schub h k , Y das gefilterte Detektorbild, X das ungefilterte 
Detektorbild, R 7 das Rampenfilter und M k die Maskierung je- 

weils bezogen auf die Indexzahl k des Liegenvorschubes h k 
darstellen . 

Weiterhin wird als besonders vorteilhafte Ausflihrung vorge- 
schlagen, eine Glattung beziiglich der Maskierung vorzunehmen. 
Bevorzugt kann hier die Maskierungsoperation so gestaltet 



\ 
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werden, dass es im Ubergangsbereich eine sin 2 -Funktion ent- 
spricht. Auf jecien Fall muss fur eine Glattung eine stetige 
Funktion verwendet werden. Die Glattung kann dabei eindimen- 
sional in eine Richtung oder zweimal hintereinander eindimen- 
5 sional in zwei Richtungen erfolgen, wobei diese Richtungen 
vorzugsweise senkrecht zueinander stehen. Eindimensionale 
Glattungen bergen dabei den Vorteil wenig rechenintesiv zu 
sein und damit schnell durchgefuhrt werden zu konnen. Alter- 
nativ kann auch eine zweidimensionale Glattung durchgefuhrt 
10 werden, die zwar bezuglich Rechenleistung aufwandiger ist, 

sich jedoch bei einetn Durchgang urn einiges effektiver gestal- 




tet . 



Neben dem oben geschilderten Verfahren schlagen die Erfinder 
• 15 auch ein CT-Gerat zur Abtastung eines Untersuchungsob j ekts 

vor, welches mit einem von mindestens einem Fokus ausgehenden 
Strahlenbundel und mit einem flachig ausgebildeten Detektor- 
array mit einer Vielzahl von verteilten Detektorelementen zum 
Detektieren der Strahlen des Strahlenbundels ausgestattet 
2 0 ist, wobei der mindestens eine Fokus relativ zu dem Untersu- 
chungsobjekt auf mindestens einer das Untersuchungsob jekt um- 
laufenden Fokusbahn mit gegenuberliegendem Detekt or array be- 
wegt. Aulierdem soil das CT-Gerat zumindest Mittel zur Samm- 
lung von Detektordaten, Filterung und Ruckpro j ektion aufwei- 
sen, wobei die Mittel zur Filterung derart gestaltet sind, 
r ' 7 dass das oben beschriebene Verfahren durchgefuhrt wird. Vor- 
zugsweise werden die genannten Mittel zur Filterung zumindest 
teilweise durch Programme oder Programm-Module" verwirklicht 
werden . 

30 

Im folgenden wird die Erfindung anhand eines/der bevorzugten 
Ausf iihrungsbeispiele (s) mit Hilfe der Figuren/der Figur naher 
beschrieben. Es zeigen im einzelnen: 
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Figur 1 : 
Figur 2 : 



Mehrzeilen-CT in perspektivischer Darstellung der 
Abtasteinheit und Peripheriegeraten; 
Langsschnitt durch das Mehrzeilen-CT aus Figur 1; 



\ 
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Figur 4 : 
Figur 5 : 



Figur 3 : 



Darstellung der Au f nahmegeome trie eines Spiral - 
CT-Gerates; 

Veranschauli chung einer Maskierungsoperat ion ; 
Beispiel von ungef ilterten und gefilterten Detek- 
tordaten anhand einer Computersimulat ion . 



10 
15 



20 
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Die Figuren 1 und 2 zeigen eine teils perspektivische Dar- 
stellung eines zur Durchfuhrung des erf indungsgemaften Verfah- 
rens geeigneten Mehrzeilen-CT-Gerates der dritten Generation. 
Die mit 1 bezeichnete Messanordnung (=Gantry) weist eine 
Rontgenstrahlenquelle 2 mit einer dieser vorgelagerten quel- 
lennahen Strahlenblende 3 auf. Das flachig ausgebildete Array 
mit mehreren Zeilen und Spalten von Detektorelementen bilde- 
tet das Detektorsystem 5 und ist mit einer diesem vorgelager- 
ten detektornahen Strahlenblende 6 in der Figur 2 im Schnitt 
dargestellt. In der Figur 1 sind zur besseren Ubersicht lich- 
keit nur acht Zeilen L x bis L Q von Detektorelementen 4 darge- 
stellt. Das Detektorsystem 5 kann jedoch auch eine andere, v 
vorzugsweise grofiere Zeilenzahl aufweisen, ohne den Rahmen 
der Erfindung zu verlassen. Ebenso ist auch eine andere fla- 
chige Anordnung der Detektoren moglich. 

Die Rontgenstrahlenquelle 2 mit der Strahlenblende 3 einer- 
seits und das Detektorsystem 5 mit der Strahlenblende 6 ande- 
rerseits sind an einem Drehrahmen 7 einander derart gegenu- 
berliegend angebracht, dass ein im Betrieb des CT-Gerats von 
der Rontgenstrahlenquelle 2 ausgehendes, durch die einstell- 
bare Strahlenblende 3 eingeblendetes, pyramidenf ormiges Ront- 
genstrahlenbundel, dessen Randstrahlen mit 8 bezeichnet sind, 
auf das Detektorsystem 5 auftrifft. Dabei ist die Strahlen- 
blende 6 dem mittels der Strahlenblende 3 eingestellten Quer- 
schnitt des Rontgenstrahlenbundels entsprechend so einge- 
stellt, dass entsprechend unterschiedlicher Betriebsmodi nur 
derjenige Bereich des Detektorsystems 5 freigegeben ist, der 
von dem Rontgenstrahlenbundel unmittelbar getroffen wird. In 
den Figuren 1 und 2 werden nur acht Zeilen von Detektorele- 
menten 4, verwendet, die weiteren punktiert angedeuteten Zei- 
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len sind von der Strahlenblende 6 abgedeckt und daher nicht 
aktiv. 

Das Rontgenstrahlenbundel weist einen Conewinkel <p auf , bei 
dem es sich urn den Of f nungswinkel des Rontgenstrahlenbundels 
in einer die Systemachse Z und den Fokus F enthaltenden Ebene 
handelt. Der Of f nungswinkel des Rontgenstrahlenbundels in ei- 
ner rechtwinklig zu der Systemachse Z liegenden und den Fokus 
F enthaltenden Ebene ( Facheroff nungswinkel ) betragt 

Der Drehrahmen 7 kann mittels einer Antriebseinrichtung 22 um 
eine mit Z bezeichnete Systemachse in Rotation versetzt wer- 
den. Die Systemachse Z verlauft parallel zu der z-Achse eines 
in der Figur 1 dargestellten und raumlich rechtwinkligen Ko- 
ordinatensystems . 

< 

Die Spalten des Detektorsystems 5 verlaufen ebenfalls in 
Richtung der z-Achse, wahrend die Zeilen quer zu der System- 
achse Z verlaufen. 

Um ein Untersuchungsobj ekt , z.B. einen Patienten, in den 
Strahlengang des Rontgenstrahlenbundels bringen zu konnen, 
ist eine Lagerungsvorrichtung 9 vorgesehen, die parallel zu 
der Systemachse Z, also in Richtung der z-Achse, verschiebbar 
ist, und zwar derart, dass eine Synchronisation zwischen der 
- f * Rotationsbewegung des Drehrahmens 7 und der Translat ionsbewe- 

gung der Lagerungsvorrichtung in dem Sinne vorliegt. Das Ver- 
haltnis von Translations- zu Rotationsgeschwindigkeit ist 
einstellbar, indem ein gewunschter Wert fur den Vorschub h 
3 0 der Lagerungsvorrichtung 9 pro Umdrehung des Drehrahmens ge- 
wahlt wird. 

Es kann also ein Volumen eines auf der Lagerungsvorrichtung 9 
befindlichen Untersuchungsobj ekts im Zuge einer Volumenabtas- 
35 tung untersucht werden, wobei die Volumenabtastung in Form 
einer Spiralabtastung in dem Sinne vorgenommen werden kann, 
dass unter gleichzeitiger Rotation der Messeinheit 1 und 
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Translation der Lagerungsvorrichtung 9 mittels der Messein- 
heit pro Umlauf der Messeinheit 1 eine Vielzahl von Projekti- 
onen aus verschiedenen Proj ektionsrichtungen auf genommen 
wird. Bei der Spiralabtastung bewegt sich der Fokus F der 
5 Rontgenstrahlenquelle relativ zu der Lagerungsvorrichtung 9 
auf einer Spiralbahn S urn das Untersuchungsobj ekt . 

Wahrend der Spiralabtastung werden aus den Detektorelementen 
jeder aktiven Zeile des Detektorsystems 5 Messdaten parallel 
10 ausgelesenen, die einzelnen Projektionen in Kegelstrahlgeo- 

metrie entsprechen, in einer Datenauf bereitungseinheit 10 ei- 
# _Vv. ner Digital/Analog-Wandlung unterzogen, serialisiert und an 
einen Bildrechner 11 iibertragen. 

■ 15 Nach einer Vorverarbeitung der Messdaten in einer Vorverar- 
beitungseinheit 12 des Bildrechners 11 gelangt der resultie- 
rende Datenstrom zu einer Volumenbildrekonstruktionseinheit 
13, die aus den Messdaten Bilder des gewiinschten Volumens des 
Untersuchungsob j ekts durch die spater beschriebene Rechenope- 
20 ration rekonstruiert . 

Die CT-Bilder setzen sich aus gitterartig zusammengeset zten 
Voxeln (Voxel= volume element) zusammen, wobei die Voxel dem 
jeweiligen Bildvolumen zugeordnet sind, jedem Voxel eine CT- 
yjj&~> Zahl in Hounsfield Units (HU) zugeordnet ist und die einzel- 
* A nen Voxel entsprechend einer CT-Zahl/Grauwertskala in einem 

ihrer jeweiligen CT-Zahl entsprechenden Grauwert dargestellt 
werden konnen. Das Bildvolumen kann auf unterschiedliche Wei- 
se visualisiert werden. In einer einfachen Form konnen belie- 

30 bige Ebenen in beliebiger Orientierung als Schnittbilder dar- 
gestellt werden. Es gibt jedoch auch komplexere Verfahren, 
die das gesamte Volumen visualisieren . Beispielhaft sei hier 
das Shaded Surface Display (SSD) und das Volume Rendering 
(VR) genannt . Das SSD berechnet zu einem einstellbaren 

35 Schwellenwert eine Isoflache, bei der das Bildvolumen den 
Schwellenwert besitzt. Die Isoflache kann mit Methoden der 
Computergraphik auf einem Bildschirm dargestellt werden. Bei 
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dem VR werden jedem Voxel gemafi seinem Wert einstellbare op- 
tische Eigenschaf ten wie Opazitat und Farbe zugeordnet. Mit 
Methoden der Computergraphik werden kunstliche Ansichten auf 
das solchermalien definierte Objekt berechnet. 

Aus den von der Bildrekonstruktionseinheit 13 rekonstruierten 
Volumendaten werden von einer Visualisierungseinheit 15 Bil- 
der berechnet und auf einer an den Bildrechner 11 angeschlos- 
senen Anzeigeeinheit 16, z.B. einem Monitor, dargestellt. 

Die Rontgenstrahlenquelle 2, beispielsweise eine Rontgenroh- 
re, wird von einer Generatoreinheit 17 mit den notwendigen 
Spannungen und Stromen, beispielsweise der Rohrenspannung £/, 
versorgt . Um diese auf die jeweils notwendigen Werte einstel- 
len zu konnen, ist der Generatoreinheit 17 eine Steuereinheit 
18 mit Tastatur 19 zugeordnet, die die notwendigen Einstel- 
lungen gestattet. 

Auch die sonstige Bedienung und Steuerung des CT-Gerates er- 
folgt mittels der Steuereinheit 18 und der Tastatur 19, was 
dadurch veranschaulicht ist, dass die Steuereinheit 18 mit 
dem Bildrechner 11 verbunden ist. 

Unter anderem kann die Anzahl der aktiven Zeilen von Detek- 
torelementen 4 und damit die Position der Strahlenblenden 3 
und 6 eingestellt werden, wozu die Steuereinheit 18 mit den 
Strahlenblenden 3 und 6 zugeordneten Verstelleinheiten 20 und 
21 verbunden ist. Weiter kann die Rotat ionszeit r eingestellt 
werden, die der Drehrahmen 7 fur eine vollstandige Umdrehung 
benotigt, was dadurch veranschaulicht ist, dass die dem Dreh- 
rahmen 7 zugeordnete Antriebseinheit 22 mit der Steuereinheit 
18 verbunden ist . 

In der Figur 3 ist die Auf nahmegeometrie des Spiral -CT- 
Gerates aus den Figuren 1 und 2 schematisch dargestellt, wo- 
bei hier - sowie in der Figur 4 - der Einfachheit halber der 
Detektor nur planar widergegeben ist. Die Strahlenquelle be- 
wegt sich mit ihrem Fokus F entlang einer Spiralbahn S mit 
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der Steigung h, die dem Liegenvorschub entspricht, rotierend 
um ein lang gestrecktes Untersuchungsobj ekt P und in Richtung 
Systemachse Z fort. Hierbei wird auf dem gegenuberliegenden 
Detektor D die Strahlenintensitat der das Untersuchungsobj ekt 
5 P durchdringenden Strahlen gemessen und gesammelt. Innerhalb 
des Untersuchungsobj ektes P ist ein zu rekonstruierendes Vo- 
lumen V R eingebettet . 



Zur Rekonstruktion eines Volumenbildes aus den gemessenen Da- 
10 ten mussen diese gefiltert und die gefilterten Daten an- 
schlieftend zur Erzeugung eines Volumenbildes des Untersu- 
chungsobj ekts dreidimensional riickpro j iziert werden, wobei 
dann das Volumenbild aus den Daten gewonnene Absorptionswerte 
der zu dem Volumen des Untersuchungsobj ekts gehorigen Voxel 
- 15 fur -die Strahlung des Strahlenbundels reprasentiert . 

Er f indungsgemafi kann die Filterung durch mehrfache Anwendung 
eines Rampenf liters R r und einer Maskierungsoperat ion M auf 

ein Proj ektionsbild in unterschiedlicher Reihenfolge erfol- 
2 0 gen. 

Die Figur 4 veranschaulicht eine' solche Maskierungsoperat ion 
liber die Detektordaten . Hintergrund ist, dass bei der Messung 
nur bestimmte Daten des messenden Detektors D zur Rekonstruk- 
v^2 A 5 tion von Interesse sind und zusatzliche Inf ormationen storend 



wirken wurden. Die Daten uber den Detektor D werden also. der- 
art durch eine Maskierungsoperat ion M maskiert, dass aus- 
schlieftlich die innerhalb der Maskierungsgrenzlinien V Q ben und 
V U nten liegenden Daten erhalten bleiben. Hierbei werden die 
3 0 obere und untere Grenze der Maskierung durch die Proj ekt ion 
der Spiralbahn S ausgehend vom Fokus F auf den Detektor ge- 
bildet. Zusatzlich werden die gemessenen Detektordaten auch 
noch einer zweidimensionalen Rampenf ilterung in Richtung der 
auf den Detektor projizierten Spiraltangente t unterzogen. 

35 



Erf indungsgemali wird dann vor der Rekonstruktion der Schnitt- 
bilder die Filterung der Detektordaten beispielsweise durch 
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die zweifache und umgekehrte Anwendung einer Rampenf ilter- 
und Maskierungsoperation auf die ungef ilterten Detektordaten 
durchgef uhrt . Eine Normierung entsprechend der Anzahl der O- 
perationen ist selbstverstandlich . Hierdurch ist eine sehr 
5 einfache, jedoch der exakten rechnerischen Losung sehr nahe 
kommende , Fi 1 1 ervor schri f t gef unden . 

Zur Veranschaulichung ist zusatzlich die erf indungsgemafie 
Filteroperation in der Figur 5 dargestellt. Diese zeigt links 
10 zunachst ein ungef iltertes Detektorbild X, auf welches - mit 

Pfeil nach oben - die Maskierungsoperation M gefolgt vom 
\ Rampenf ilter R t angewendet wurde, wodurch sich das Detektor- 
fPfy^ bild R t MX ergibt . Mit Pfeil nach unten ist auch das Ergebnis 
der gleichen, jedoch in umgekehrter Reihenfolge ausgefiihrten 
. 15 Operationen dargestellt, die das Detektorbild MRfX ergeben. 

Anschlieftend werden beide Ergebnisse addiert und normiert, so 
dass sich als Endergebnis durch einfache Rechenoperat ionen 
das fertige, gefilterte Detektorbild Y ergibt. 

20 Insgesamt wird also beim erf indungsgemaften Verfahren die be- 
kannte Filterung der Detektordaten bei an sich bekannten Re- 

konstrukt ionsverf ahren von Volumenbildnern durch mehrfache 
Anwendung eines Rampenf ilters R 7 und einer Maskierungsopera- 
tion M auf ein Proj ektionsbild in unterschiedlicher Reihen- 
i 'il 2 * 5 folge ersetzt und damit ein Verfahren zur Bilderstellung fur 
m ~*jr einen mehrzeiligen Spiral -Computertomographen offenbart, das 
. einerseits die guten Eigenschaf ten der exakten Verfahren be- 
zuglich der Bildqualitat beibehalt, andererseits jedoch durch 
die Verwendung einer Naherung den Berechnungsauf wand dras- 
30 tisch verringert und damit auch die Flexibilitat der Berech- 
nung erhoht . 
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Patentanspruche 

Verfahren zur Bilderstellung fur einen mehrzeiligen 
Spiral - Comput e r t omogr aphen , aufweisend die Verfahrens- 
schritte : 

ein Untersuchungsobj ekt wird durch ein rotierendes und 
sich in Richtung der Rotationsachse bewegendes Strah- 
lenbundel bezuglich seines Absorptionsverhaltens abge- 
tastet und die gemessenen Absorpt ionsdaten werden ge- 
sammelt , 

zur Rekonstruktion eines Volumenbildes aus den gemesse- 
nen Daten werden diese gefiltert und die gefilterten 
Daten anschlieftend zur Erzeugung eines Volumenbildes 
des Untersuchungsob j ekts dreidimensional ruckpro j i- 
ziert, wobei das Volumenbild aus den Daten gewonnene 
Absorptionswerte der zu dem Volumen des Untersuchungs- 
objekts gehorigen Voxel fur die Strahlung des Strahlen- 
bundels reprasentiert , 

dadurch gekennzeichnet, dass 
die Filterung durch mehrfache Anwendung eines Rampen- 
f ilters R 7 und einer Maskierungsoperation M auf ein 

Projektionsbild in unterschiedlicher Reihenfolge er- 
folgt . 

Verfahren gemaft dem voranstehenden Patentanspruch 1, 
dadurch gekennzeichnet, dass das Rampenfilter R 7 

und die Maskierungsoperation M abwechselnd angewendet 
wird . 

30 3. Verfahren gemaft einem der voranstehenden Patentanspru- 
che 1 bis 2, dadurch gekennzeichnet, dass die 
Filterung mit der Vorschrift: 

Y = ^(R 7 M + MR 7 )X 

durchgefuhrt wird, worin X das ungef ilterte Detektor- 
35 bild, R 7 die Rampenf ilterung in Richtung der Projektion 

der Spiraltangente und M die Maskierung darstellen. 
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4. Verfahren gemafi einem der voranstehenden Patentanspru- 
che 1 bis 2, dadurch gekennzeichnet , dass ein 
Liegenvorschub h k verwendet wird, dessen Wert der Be- 
dingung 



gentigt, wobei k eine natiirliche Zahl mit k > 0 ist und 
h m ax dem Liegenvorschub entspricht, bei dem eine durch 
die Maskierung ausgezeichnete Flache gerade noch den 
Detektor abdeckt . 

5. Verfahren gemaft dem voranstehenden Patentanspruch 4, 
dadurch gekennzeichnet, dass die Filterung mit 
der Vorschrift: 

Y = — l —(R 7 M. +M,R 7 )X 

durchgefilhrt wird, mit k der Indexzahl aus dem Liegen- 
vorschub hk, Y dem gef ilterten Detektorbild, X dem un- 
gefilterten Detektorbild, R 7 dem Rampenfilter und M* 

der Maskierung. 

6. Verfahren gemaB einem der voranstehenden Patentanspru- 
che 1 bis 5, dadurch gekennzeichnet, dass die 
Maskierung zumindest in ihrem Ubergangsbereich durch 
eine stetige Funktion geglattet wird. 

7. Verfahren gemaB dem voranstehenden Patentanspruch 6, 
dadurch gekennzeichnet, dass die Glattung der 
ursprunglich unstetigen Maskierung eindimensional er- 
folgt 

8. Verfahren gemafi dem voranstehenden Patentanspruch 7, 
dadurch gekennzeichnet, dass die eindimensionale 
Glattung mehrfach in unterschiedliche, vorzugsweise 
senkrecht zueinander verlaufende, Richtungen erf olgt . 

9. Verfahren gemali dem voranstehenden Patentanspruch 6, 
dadurch gekennzeichnet, dass die Glattung der 



\ 
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ursprunglich unstetigen Maskierung zweidimensional er- 
folgt . 

CT-Gerat zur Abtastung eines Untersuchungsob j ekts mit 
einem von mindestens einem Fokus ausgehenden Strahlen- 
biindel und mit einem flachig ausgebildeten Detektorar- 
ray mit einer Vielzahl von verteilten Detektorelementen 
zum Detektieren der Strahlen des Strahlenbundels , wobei 
der mindestens eine Fokus relativ zu dem Untersuchungs- 
objekt auf mindestens einer das Untersuchungsob j ekt um- 
laufenden Fokusbahn mit gegeniiberliegendem Detektorar- 
ray bewegt, wobei zumindest Mittel zur Sammlung von De- 
tektordaten, Filterung und 3D Riickpro j ektion vorgesehen 
sind und die Mittel zur Filterung derart gestaltet 
sind, dass das Verfahren gemafi einem der Anspruche 1 
bis 9 durchgefuhrt wird. 

11. CT-Gerat gemafi Anspruch 10, dadurch gekenn- 

zeichnet, dass die genannten Mittel zur Filterung 
2 0 zumindest teilweise durch Programme oder Programm- 

Module verwirklicht werden. 
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Zusammenf as sung 

Verfahren zur Bilderstellung fur einen mehrzeiligen Spiral- 
Computertomographen und Computer-Tomographie-Gerat zur Durch- 
fuhrung dieses Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bilderstellung fur 
einen mehrzeiligen Spiral -Computertomographen, wobei ein Un- 
tersuchungsobj ekt wird durch ein rotierendes und sich in 
Richtung der Rotationsachse bewegendes Strahlenbundel beziig- 
lich seines Absorpt ionsverhaltens abgetastet und die gemesse- 
nen Absorpt ionsdat en werden gesammelt, und zur Rekonstruktion 
von Volumenbildern aus den gemessenen Daten werden diese ge- 
filtert und die gefilterten Daten anschlieBend zur Erzeugung 
wenigstens eines Volumenbildes des Unt er sue hung sobj ekts drei- 
dimensional riickproj iziert , wobei das Volumenbild aus den Da- 
ten gewonnene Absorptionswerte der zu der Schicht des Unter- 
suchungsob j ekts gehorigen Voxel fur die Strahlung des Strah- 
lenbundels reprasentiert . Aufterdem betrifft die Erfindung ein 
CT-Gerat zur Durchfuhrung dieses Verfahrens mit entsprechen- 
den Mitteln zur Bilderstellung. 

Beim erfindungsgemaBen Verfahren erfolgt die Filterung durch 
mehrfache Anwendung eines Rampenf ilters R 7 und einer Maskie- 

rungsoperation M auf ein Projektionsbild in unterschiedli - 
cher Reihenfolge. Das erf indungsgemafie CT-Gerat weist Mittel 
auf, die zur Filterung mehrfach einen Rampenf ilter und eine 
Maskierungsoperation auf ein Projektionsbild in unterschied- 
licher Reihenf olge anwenden . 



Figur 5 
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